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Abstract. We study the densities of uniform random walks in the plane. A spe- 
cial focus is on the case of short walks with three or four steps and less completely 
those with five steps. As one of the main results, we obtain a hypergeometric rep- 
resentation of the density for four steps, which complements the classical elliptic 
representation in the case of three steps. It appears unrealistic to expect similar 
results for more than five steps. New results are also presented concerning the 
moments of uniform random walks and, in particular, their derivatives. Relations 
with Mahler measures are discussed. 



An n-step uniform random walk is a walk in the plane that starts at the origin 
and consists of n steps of length 1 each taken into a uniformly random direction. 
The study of such walks largely originated with Pearson more than a century ago 
[PeaOSb, PeaOSa, Pea06] who posed the problem of determining the distribution of 
the distance from the origin after a certain number of steps. In this paper, we study 
the (radial) densities Pn of the distance travelled in n steps. This continues research 
commenced in [BNSW09, BSWll] where the focus was on the moments of these 
distributions: 



The densities for walks of up to 8 steps are depicted in Figure 1. As established 
by Lord Rayleigh [Ray05], pn quickly approaches the probability density ^e""^ for 
large n. This limiting density is superimposed in Figure 1 for n ^ 5. 



Closed forms were only known in the cases n = 2 and n = 3. The evaluation, for 

^ X ^ 2, 



with an appendix by: DON ZAGIER 



1. Introduction 




(1) 



P2{x) 



2 




2010 Mathematics Subject Classification. Primary 60G50; Secondary 33C20, 34M25, 44A10. 



2 



JONATHAN M. BORWEIN, ARMIN STRAUB, JAMES WAN, AND WADIM ZUDILIN 




is elementary. On the other hand, the density ^3(0:) for ^ x ^ 3 can be expressed 
in terms of elhptic integrals by 



(2) 




X + 1)3(3 -X 

16x 




see, e.g., [Pea06]. One of the main results of this paper is a closed form evaluation 
of p4 as a hypergeometric function given in Theorem 8. In (20) we also provide a 
single hypergeometric closed form for ^3 which, in contrast to (2), is real and valid 
on all of [0, 3]. For convenience, we list these two closed forms here: 



(3) 



(4) 



P4(x) 




We note that while Maple handled these well to high precision, Mathematica strug- 
gled, especially with the analytic continuation of the 3F2 when the argument is 
greater than 1. 
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A striking feature of the 3- and 4-step random walk densities is their modularity. 
It is this circumstance which not only allows us to express them via hypergeometric 
series, but also makes them a remarkable object of mathematical study. 

This paper is structured as follows: In Section 2 we give general results for the 
densities Pn and prove for instance that they satisfy certain linear differential equa- 
tions. In Sections 3, 4, and 5 we provide special results for ps, p4, and ps respectively. 
Particular interest is taken in the behaviour near the points where the densities fail 
to be smooth. In Section 6 we study the derivatives of the moment function and 
make a connection to multidimensional Mahler measures. Finally in Section 7 we 
provide some related new evaluations of moments and so resolve a central case of an 
earlier conjecture on convolutions of moments in [BSWll]. 

The amazing story of the appearance of Theorem 4 is worth mentioning here. The 
theorem was a conjecture in an earlier version of this manuscript, and one of the 
present authors communicated it to D. Zagier. That author was surprised to learn 
that Zagier had already been asked for a proof of exactly the same identities a little 
earlier, by P. Djakov and B. Mityagin. 

Those authors had in fact proved the theorem already in 2004 (see [DM04, The- 
orem 4.1] and [DM07, Theorem 8]) during their study of the asymptotics of the 
spectral gaps of a Schrodinger operator with a two-term potential — their proof was 
indirect, so that we should never have come across the identities without the acci- 
dent of asking the same person the same question! Djakov and Mityagin asked Zagier 
about the possibility of a direct proof of their identities (the subject of Theorem 4), 
and he gave a very neat and purely combinatorial answer. It is this proof which is 
herein presented in the Appendix. 

We close this introduction with a historical remark illustrating the fascination 
arising from these densities and their curious geometric features. H. Pettis devotes 
the entire paper [Pet63] to proving that ps is not linear on the initial interval [0, 1] 
as ruminated upon by Pearson [Pea06]. This will be explained in Section 5. 

2. The densities pn 

It is a classical result of Kluyver [Klu06] that Pn has the following Bessel integral 
representation: 



Here J,y is the Bessel function of the first kind of order u. 

Remark 1. Equation (5) naturally generalizes to the case of nonuniform step lengths. 
In particular, for n = 2 and step lengths a and b we record (see [Wat41, p. 411] or 



(5) 
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[Hug95, 2.3.2]; the result is attributed to Sonine) that the corresponding density is 

/■oo 

P2{x;a,b)= / xtJo{xt)Jo{at)Jo{bt)dt 
Jo 

(6) = , ^ = 

7ry/{{a + by-x^){x^-{a-by) 

foT \a — b\ ^ X ^ a + b and p2{x; a,b) = otherwise. Observe how (6) speciahzes to 
(1) in the case a = b = 1. 

In the case n = 3 the density psix; a, b, c) has been evaluated by Nicholson [Wat41, 
p. 414] in terms of elliptic integrals directly generalizing (2). The corresponding 
extensions for four and more variables appear much less accessible. 

It is visually clear from the graphs in Figure 1 that pn is getting smoother for 
increasing n. This can be made precise from (5) using the asymptotic formula for Jq 
for large arguments and dominated convergence: 

Theorem 1. For each integer n ^ 0, the density Pn+4, is \n/2\ times continuously 
differentiable. 

On the other hand, we note from Figure 1 that the only points preventing p„ from 
being smooth appear to be integers. This will be made precise in Theorem 2. 

To this end, we recall a few things about the s-th moments Wn{s) of the density 
Pn which are given by 



(7) Wnis) = / X'pnix)dx 



[0,1]- 



k=l 



e 



da;. 



Starting with the right-hand side, these moments had been investigated in [BNSW09, 
BSWll]. There it was shown that Wn{s) admits an analytic continuation to all of 
the complex plane with poles of at most order two at certain negative integers. In 
particular, 1^3(5) has simple poles at s = —2, —4, —6, . . . and ^4(5) has double poles 
at these integers [BNSW09, Thm. 6, Ex. 2 & 3]. 
Moreover, from the combinatorial evaluation 



Wn{2k) = V , 

aiH han = '" 



k 



for integers ^ it followed that Wn{s) satisfies a functional equation, as in 
[BNSW09, Ex. 1], coming from the inevitable recursion that exists for the right- 
hand side of (8) . For instance, 

(s + AfW^{s + 4) - 2(55^ + 30s + AQ)W^{s + 2) + 9(s + 2fW^{s) = 0, 

and equation (9) below. 
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The first part of equation (7) can be replirased as saying tliat Wn{s—1) is tlie Mellin 
transform of p„ ([ML86]). We denote this by Wn{s — 1) = s]. Conversely, the 

density Pn is the inverse Mellin transform of Wn{s — 1). We intend to exploit this 
relation as detailed for n = 4 in the following example. 

Example 1 (Mellin transforms). For n = 4, the moments ^4(5) satisfy the func- 
tional equation 

(9) {s + AfWiis + 4) - 4(s + 3)(5s2 + 30s + 48)^4(5 + 2) + 64(s + 2fW4{s) = 0. 

Recall the following rules for the Mellin transform: if F{s) = Ai [/; s\ then in the 
appropriate strips of convergence 

• M [a;'^/(a;);s] =F(s + /i), 

• M[DJ{x)-s\ = -{s-l)F{s-l). 

Here, and below, denotes differentiation with respect to x, and, for the second 
rule to be true, we have to assume, for instance, that / is continuously different iable. 

Thus, purely formally, we can translate the functional equation (9) of W4 into the 
differential equation ^4 ■ Pi{x) = where ^4 is the operator 

(10) A4 = x\e + I f - Ax^eihe^ + 3) + 64(^ - 1)^ 

(11) ={x-A){x- 2)x^{x + 2)(x + A)Dl + Qx^ [x"^ - lO) Dl 

+ X (7x^ - 32x2 + 64) D.^ + (x^ - S) (x^ + 8) . 

Here 9 = xD^. However, it should be noted that is not continuously differentiable. 
Moreover, p4(x) is approximated by a constant multiple of a/4 — x as x — 4~ (see 
Theorem 5) so that the second derivative of p^ is not even locally integrable. In 
particular, it does not have a Mellin transform in the classical sense. 

Theorem 2. Let an integer n ^ 1 be given. 

• The density Pn satisfies a differential equation of order n — 1. 

• If n is even (respectively odd) then Pn is real analytic except at and the even 
(respectively odd) integers m <n. 

Proof. As illustrated for p^ in Example 1, we formally use the Mellin transform 
method to translate the functional equation of Wn into a differential equation An ■ 
y{x) = 0. Since pn is locally integrable and compactly supported, it has a Mellin 
transform in the distributional sense as detailed for instance in [ML86]. It follows 
rigorously that Pn solves A„ ■ y{x) = in a distributional sense. In other words, p„ 
is a weak solution of this differential equation. The degree of this equation is n — 1 
because the functional equation satisfied by Wn has coefficients of degree n — 1 as 
shown in [BNSW09, Thm. 1]. 
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The leading coefficient of the differential equation (in terms of as in (11)) turns 
out to be 

(12) x"-^ JJ {x^-m") 

2\{m-n) 

where the product is over the even or odd integers 1 < m < n depending on whether 
n is even or odd. This is discussed below in Section 2.1. 

Thus the leading coefficient of the differential equation is nonzero on [0, n] except 
for and the even or odd integers already mentioned. On each interval not containing 
these points it follows, as described for instance in [Hor89, Cor. 3.1.6], that p„ is 
in fact a classical solution of the differential equation. Moreover the analyticity of 
the coefficients, which are polynomials in our case, implies that p„ is piecewise real 
analytic as claimed. □ 

Remark 2. It is one of the basic properties of the Mellin transform, see for instance 
[FS09, Appendix B.7], that the asymptotic behaviour of a function at zero is deter- 
mined by the poles of its Mellin transform which lie to the left of the fundamental 
strip. It is shown in [BNSW09] that the poles of Wn{s) occur at specific negative 
integers and are at most of second order. This translates into the fact that p„ has an 
expansion at as a power series with additional logarithmic terms in the presence 
of double poles. This is made explicit in the case of p4 in Example 4. 

2.1. An explicit recursion. We close this section by providing details for the claim 
made in (12). Recall that the even moments fn{k) '■= Wn{2k) satisfy a recurrence of 
order A := [?t-/2] with polynomial coefficients of degree n — 1 (see [BNSW09]). An 
entirely explicit formula for this recurrence is given in [Ver04]: 



E n(-"<)(«+i-".)(j^)' 

i=l 



fnik-j) = 



Theorem 3. 

(13) E 

where the sum is over all sequences ai, . . . ,aj such that ^ ai ^ n and aj+i ^ ctj — 2. 

Observe that (12) is easily checked for each fixed n by applying Theorem 3. We 
explicitly checked the cases n ^ 1000 (using a recursive formulation of Theorem 3 
from [Ver04]) while only using this statement for n ^ 5 in this paper. The fact that 
(12) is true in general is recorded and made more explicit in Theorem 4 below. 

For fixed n, write the recurrence for fn{k) in the form ^^Z^k^qj^K) where qj 
are polynomials and K is the shift k — )■ k + 1. Then is the characteristic 
polynomial of this recurrence, and, by the rules outlined in Example 1, we find that 
the differential equation satisfied by Pn{x) is of the form g„_i(x^)0"~^ + ■ ■ • ? where 
9 = xD^ and the dots indicate terms of lower order in 6. 
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We claim that the characteristic polynomial of the recurrence (13) satisfied by 

fn{k) is n 



2\{m-n) 



[x — m ) where the product is over the integers 1 < m < n such 
that m = n modulo 2. This implies (12). By Theorem 3 the characteristic polynomial 
is 



(14) 



A 

E 

i=o 



a,; 



ai,...,aj i=l 



X 



A-i 



where A = \n/2] and the sum is again over all sequences ai, . . . ,aj such that ^ 
ai ^ n and Oj+i ^ — 2. The claimed evaluation is thus equivalent to the following 
identity, first proven by P. Djakov and B. Mityagin [DM04, DM07]. Zagier's more 
direct and purely combinatorial proof is given in the Appendix. 

Theorem 4. For all integers n,j ^ 1, 



(15) 



0^mi,...,'mj<n/2 1=1 



l^Ja^ . . i=\ 



3. The density 



The elliptic integral evaluation (2) of ps is very suitable to extract information 
about the features of ps exposed in Figure 1(a). It follows, for instance, that has 
a singularity at 1. Moreover, using the known asymptotics for K{x), we may deduce 
that the singularity is of the form 



(16) 



P3(.) = ^l0g (-^ 



1 



+ 0(1) 



as X — 7- 1. 

We also recall from [BSWll, Ex. 5] that has the expansion, valid for ^ x ^ 1, 

2r /tn2A; 



(17) 
where 
(18) 



7T 



k=0 



Wsi2k) = J2 



j=Q 



V3 

j 



is the sum of squares of trinomials. Moreover, we have from [BSWll, Eqn. 29] the 
functional relation 



(19) 



Ax 



(3-x)(x + l) 



P3 



X 



X 
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SO that (17) determines completely and also makes apparent the behaviour at 3. 
We close this section with two more alternative expressions for p^. 

Example 2 (Hypergeometric form for p^). Using the techniques in [CZIO] we can 
deduce from (17) that 



2 V3x ^ fl 2 ^ x2 (9 



X 



2\2 



which is found in a similar but simpler way than the hypergeometric form of p4 given 
in Theorem 8. Once obtained, this identity is easily proven using the differential 
equation from Theorem 2 satisfied by p^. From (20) we see, for example, that 

PsiVsr = i,w,{-i). 

Example 3 (Iterative form for p^). The expression (20) can be interpreted in terms 
of the cubic AGM, AG3, see [BB91], as follows. Recall that AG3(a,6) is the limit of 
iterating 

_ a„ + 2bn 

O'n+l — ^ ; On+l 

beginning with qq = a and 60 = b. The iterations converge cubically, thus allowing 
for very efficient high-precision evaluation. On the other hand, 

1 fl 2 3 

AG^-^^^U'3''''~' 

so that in consequence of (20), for ^ x ^ 3, 

(21) mix) 2^ 




AG3(3 + x2,3|l -x2|2/2)' 

Note that P3(3) = |^ is a direct consequence of the final formula. 

Finally we remark that the cubic AGM also makes an appearance in the case 
n = 4. We just mention that the modular properties of p4 recorded in Remark 7 can 
be stated in terms of the theta functions 

(22) 0(r) = c(r) = 3 

r7(3r) r/(r) 

where rj is the Dedekind eta function defined in (42). For more information and 
proper definitions of the functions 6, c as well as a, which is related by = h^ + c^, 
we refer to [BBG94]. Ultimately we are hopeful that, in search for an analogue 
of (19) for p4, this may lead to an algebraic relation between algebraically related 
arguments of ^4. 
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4. The density 

The densities p„ are recursively related. As in [Hug95], setting = p„(x)/(27rx), 

we have that for integers n ^ 2 

(23) 4>n{x) = — / 0n-i [y/x'^ — 2x cos a + 1 ) da. 

27r Jo ^ / 

We use this recursive relation to get some quantitative information about the 
behaviour of p4 at x = 4. 



Theorem 5. As x 4 , 

^v/4^- 

TT^ 167r^ ^ ' 5127r 



P4(x) = ^V^x - ^(4 - xf'^ + 1^(4 - xfl- + O ((4 - xf'^) 



Proof. Set y = y/x"^ — 2x cos a + 1. For 2 < x < 4, 

1 /•'^ 1 /•arccos(4i^) 

04(a;) = - / 03(2/) da = - (p^^y) da 

Jo Jo 

since (^3 is only supported on [0,3]. Note that (psly) is continuous and bounded in 
the domain of integration. By the Leibniz integral rule, we can thus differentiate 
under the integral sign to obtain 

.// N 1 (x2 + 8)03(3) 1 (4f^) ,^^03(2/)^ 

(24) ^'Jx) = , i-r-iy j, + - / (x - cos(a)) da. 

vrxV(16-x2)(x2-4) vrio 2/ 

This shows that ^4, and hence P4, have a singularity at a; = 4. More specifically, 
d'Ax) = , - + 0(1) asx-^4". 

Here, we used that 03(3) = It follows that 

which, upon integration, is the claim to first order. Differentiating (24) twice more 
proves the claim. □ 

Remark 3. The situation for the singularity at a; = 2"^ is more complicated since in 
(24) both the integral (via the logarithmic singularity of 03 at 1, see (16)) and the 
boundary term contribute. Numerically, we find, as a; — 2+, 

2 
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On the other hand, the derivative of at 2 from the left is given by 



73 



TT 



3-^2 



112 
2' 3' 3 

1,1 



These observations can be proven in hindsight from Theorem 7. 







We now turn to the behaviour of at zero which we derive from the pole structure 
of W4 as described in Remark 2. 



J 


1 

I 3 
\ ^ 




-6 -4 -2 

-1 
-2 

-3 


2 -6 




(a) W4 (b) W5 

Figure 2. and W5 analytically continued to the real line. 



Example 4. From [BSWll], we know that W4 has a pole of order 2 at —2 as 
illustrated in Figure 2(a). More specifically, results in Section 6 give 



as s — )• — 2. It therefore follows that 

. . 3 



P4,[X) 



27r2 



9 
27r~2 



X log(x) + — ^ log(2)x + O(x^) 



as X — )■ 0. 







More generally, W4 has poles of order 2 at — 2A; for k a positive integer. Define s^^k 
and r4^fc by 

■S4,fc-1 , '"4,/c-l 



(25) 



W^4(S) 



+ 



+ 0(1) 



{s + 2ky s + 2k 
as s — )■ — 2A;. We thus obtain that, as x — )■ 0+, 

K-l 

P4{x) = x^^+^ (r4,fe - S4,fc log(x)) + O 



k=0 
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In fact, knowing that solves the hnear Fuchsian differential equation (10) with a 
regular singularity at we may conclude: 

Theorem 6. For small values x > 0, 

oo 

(26) p4{x) = ^ (r4,fc - S4,fc log(x)) X 

Note that 



2k+l 



k=0 



3 W4{2k) 



27r2 82^= 

as the two sequences satisfy the same recurrence and initial conditions. The numbers 
W4{2k) are also known as the Domb numbers ([BBBG08]), and their generating 
function in hypergeometric form is given in [Rog09] and has been further studied in 
[CZIO]. We thus have 



(27) 



A;=0 



SAkX 



6x 



7r2 (4 



x^ 



3-^2 



112 
3' 2' 3 

1,1 



108x2 



(x2 - 4)' 



which is readily verified to be an analytic solution to the differential equation satisfied 
by p4. 

Remark 4. For future use, we note that (27) can also be written as 



{2i 



E 

fc=0 



2k+l 



24x 



7r2 (16 



x^ 



3-^2 



112 
3' 2' 3 

1,1 



108x^ 



which follows from the transformation 

3' 2' 3 108^/ 

1,1 



(29) (1 - 4x)3F2 
given in [CZIO, (3.1)]. 



'1 - 16x) 



{16-x^f 
(1 - 16x)3F2 



112 

3' 2' 3 

1,1 



108x2 



;i -4x)= 







On the other hand, as a consequence of (25) and the functional equation (9) 
satisfied by W4, the residues r4fc can be obtained from the recurrence relation 



(30) 



128A;V4,fc = 4(2A; - l){5k^ - 5k + 2)r4,fc-i - 2{k - l)%fc_2 
+ 3 (64A;2s4,fc - (20P - 20A; + 6)54,^-1 + (k - 



with r4 _i = and r^^o = ^ log(2). 

Remark 5. In fact, before realizing the connection between the Mellin transform 
and the behaviour of at 0, we empirically found that p^ on (0, 2) should be of the 
form r(x) — s(x)log(x) where a and r are odd and analytic. We then numerically 
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determined the coefficients and observed the relation with the residues of W4 as given 



The differential equation for has a regular singularity at 0. A basis of solutions at 
can therefore be obtained via the Frobenius method, see for instance [lnc26] . Since 
the indicial equation has 1 as a triple root, the solution (27) is the unique analytic 
solution at while the other solutions have a logarithmic or double logarithmic 
singularity. The solution with a logarithmic singularity at is explicitly given in 
(34), and, from (26), it is clear that p4 on (0,2) is a linear combination of the 
analytic and the logarithmic solution. 

Moreover, the differential equation for is a symmetric square. In other words, it 
can be reduced to a second order differential equation, which after a quadratic substi- 
tution, has 4 regular singularities and is thus of Heun type. In fact, a hypergeometric 
representation of p^ with rational argument is possible. 

Theorem 7. For 2 < x < A, 



Proof. Denote the right-hand side of (31) by qi{x) and observe that the hypergeo- 
metric series converges for 2 < a; < 4. It is routine to verify that 54 is a solution of 
the differential equation ^4 ■ y{x) = given in (10) which is also satisfied by p^ as 
proven in Theorem 2. Together with the boundary conditions supplied by Theorem 



We note that Theorem 7 gives 2^16 - xV(7r2 x) as an approximation to p4{x) near 
X = 4, which is much more accurate than the elementary estimates established in 
Theorem 5. 

Corollary 1. In particular, 



in Theorem 6. 







(31) 




5 it follows that ^4 = ^4. 



□ 



(32) 




Quite marvelously, as ffist discovered numerically: 



Theorem 8. For < x < 4 



(33) 
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Proof. To obtain the analytic continuation of the 3F2 for < x < 2 we employ the 
formula [Luk69, 5.3], valid for all z, 



ai, 



. . . , ttq+l 



9+1 



r(afc) Uj^k - flfe) 



X q+lFq 



Ofc, {Ofe - 6j + l}j 

{o-k — o.j + iLyfc 



which requires the aj to not differ by integers. Therefore we apply it to 

3-^2 



2 ^ ' 2' 2 



5 7 
6' 6 



and take the limit as e — )■ 0. This ultimately produces, for z > 1, 

log(108z) 



(34) RegFa 



/I 1 1 




1 2' 2' 2 


■) 


V ^ ^ 


\ 6' 6 





+ 



2V3-Z 
1 



/ 1 1 2 




( 3' 2' 3 


V 1,1 





3-^2 



00 (1\ (1 



2V32 



n=0 



3n; 



Here if„ = 1/A; is the n-th harmonic number. Now, insert the appropriate 

argument for z and the factors so the left-hand side corresponds to the claimed 
closed form. Observing that 



{^n (2)71 (3)71 



(2n)!(3n)! 



108"(n!)2 ' 

we thus find that the right-hand side of (33) is given by — log(x)S'4(x) plus 
-2 E ''"!;fi!^'"n/!'l3„ - 2^^- - 3^^.» + 3 log(16 - x^)) 



n=0 



(16 -x2; 



where 5*4 is the solution (analytic at 0) to the differential equation for p/^ given in 
(28). This combination can now be verified to be a formal and hence actual solution 
of the differential equation for ^4. Together with the boundary conditions supplied 
by Theorem 6 this proves the claim. □ 

Remark 6. Let us indicate how the hypergeometric expression for given in The- 
orem 7 was discovered. Consider the generating series 



(35) 



k=0 
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of the Domb numbers which is just a rescaled version of (27). Corresponding to (28), 
the hypergeometric form for this series given in [Rog09] is 

108^^ 

(36) y^^,) = —^^F,— 



1-Az 



112 
3' 2' 3 

1,1 



;i - Azy 



which converges for 1^1 < 1/16. yo satisfies the differential equation ■ yo{z) 
where 



(37) 
and 6 



'dz 



Bi = 64z\9 + 1)3 - 2z{29 + 1){59^ + 59 + 2) + 9^ 
. Up to a change of variables this is (10); yo is the unique solution which 



is analytic at zero and takes the value 1 at zero; the other solutions which are not a 
multiple of yo have a single or double logarithmic singularity. Let yi be the solution 
characterized by 

(38) yi{z)-yo{z)\og{z)ezQ[[z]]. 

Note that it follows from (38) as well as Theorem 6 together with the initial values 

S4,0 

, , 3x 



and r4^o = 'S4,olog(8) that p4, for small positive argument, is given by 



(39) 



471^ 



yi 



64 ) 



If X G (2,4) and z = x^/64 then the argument t = (I'^^^yi of the hypergeometric 
function in (36) takes the values (l,oo). We therefore consider the solutions of the 
corresponding hypergeometric equation at infinity. A standard basis for these is 



(40) t-'/^F2 



111 

3' 3' 3 
2 5 
3' 6 



t-^/^3i^2 



111 
2' 2' 2 

5 1 
6' 6 



t-^/^3i^2 



2 2 2 
3' 3' 3 
I 7 
3' 6 



In fact, the second element suffices to express on the interval (2,4) as shown in 
Theorem 7. 



We close this section by showing that, remarkably, has modular structure. 

Remark 7. As shown in [CZIO] the series yo defined in (35) possesses the modular 
parameterization 

/ 77(2r)6r^(6r)6 \ _ yirMSry 
^ ' Vir^vi^ry J r/(2r)2r7(6r)2- 

Here i] is the Dedekind eta function defined as 



(42) 



^{r) = Y[{1 - g") = gV24 ^ (_i)n^„(3n+l)/2^ 
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where q = e^'^*^. Moreover, the quotient of the logarithmic solution yi defined in 
(38) and z/q is related to the modular parameter r used in (41) by 



(43) exp = 

Combining (41), (43) and (39) one obtains the modular representation 
. ^ f v(2TYr](6Tf\ 6(2r + l) , , , , , , , , 



valid when the argument of is small and positive. This is the case for r = 
—1/2+iy when y > 0. Remarkably, the argument attains the value 1 at the quadratic 
irrationality r = (-^—5/3 — l)/2 (the 5/3rd singular value of the next section). As 
a consequence, the value P4(l) has a nice evaluation which is given in Theorem 9. 



5. The density 

As shown in [BSWll], ^^5(5) has simple poles at —2, —4, . . ., compare Figure 2(b). 
We write r^^k = Res_2fc-2 for the residue of at s = —2k — 2. A surprising 
bonus is an evaluation of r^^ = ^4(1) ~ 0.3299338011, the residue at s = —2. This 
is because in general for n ^ 4, one has 

Res_2lV„+i=K+i(0)=Pn(l), 

as follows from [BSWll, Prop. 1(b)]; here p'„_,_i(0) denotes the derivative from the 
right at zero. 

Explicitly, using Theorem 8, we have, 

(45) r,o=p'5(0) = ^Re3F2(y;^ ^) . 

\ 6' 6 / 

In fact, based on the modularity of p^ discussed in Remark 7 we find: 
Theorem 9. 



(46) r5,o 



Proof. The value r = {^-5/3 - l)/2 in (44) gives the value ^4(1) = r^ Q. Applying 
the Chowla-Selberg formula [SC67, BB98] to evaluate the eta functions yields the 
claimed evaluation. □ 
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Using [BZ92, Table 4, (ii)], (46) may be simplified to 

(47) ^5,0 - — ^ 



(48) = ^ ^— ^ Kl, = ^ K,/,K,„ 



3V5 (v^- 1) ^2 _ ^ 
7r3 2 7r3 

where Ki^ and -ft's/s are the complete elliptic integral at the 15th and 5/3rd singular 
values [BB98]. 

Remarkably, these evaluations appear to extend to r^^i ^ 0.006616730259, the 
residue at s = —4. Resemblance to the tiny nome of Bologna [BBBG08] led us to 
discover — and then check to 400 places using (55) and (56) — that 

u.. , 1 13 r(^)r(^)r(A)r(A) i r(^)r(i|)r(i)r(M) 

^ ^ 1800v^ vr^ TT^ ■ 

Using (47) this evaluation can be neatly restated as 

7 13 2 1 

(50) rs,! = — r5,o 



225 ' 57r4r5,o 
We summarize our knowledge as follows: 

Theorem 10. The density is real analytic on (0, 5) except at 1 and 3 and satisfies 
the differential equation ■ p^ix) = where is the operator 

(51) As = x^(e + 1)^-^^(350^ + 42^2 ^3) 

+ x\259{e - 1)^ + 104(0 - 1)2) - (15(0 - 3)(0 - 1))' 

and 9 = xD^. Moreover, for small x > 0, 



(52) p,{x) = 
where 

{15{2k + 2){2k + 4))'r5,fc+2 = (259(2A; + 2)^ + 104(2A; + 2)^) rs.^+i 

(53) - (35(2A; + 1)^ + 42(2A; + 1)^ + 3) r^^k + {2k)%k-i 

with explicit initial values r^ _i = and r^^, r^ i given by (47) and (49) above. 

Proof. First, the differential equation (51) is computed as was that for p^, see (10). 
Next, as detailed in [BSWll, Ex. 3] the residues satisfy the recurrence relation (53) 
with the given initial values. Finally, proceeding as for (26), we deduce that (52) 
holds for small x > 0. □ 
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Numerically, the series (52) appears to converge for |x| < 3 which is in accordance 
with ^ being a root of the characteristic polynomial of the recurrence (53); see also 
(12). The series (52) is depicted in Figure 3. 



0.8 - 




0.5 1.0 1.5 2.0 2.5 



Figure 3. The series (52) (dotted) and ps. 

Since the poles of W5 are simple, no logarithmic terms are involved in (52) as 
opposed to (26). In particular, by computing a few more residues from (53), 

P5(x) = 0.329934 X + 0.00661673 + 0.000262333 + 0.0000141 185 + O(x^) 

near (with each coefficient given to six digits of precision only), explaining the 
strikingly straight shape of P5(x) on [0, 1] . This phenomenon was observed by Pearson 
[Pea06] who stated that for p5(x)/x between x = and x = 1, 

"the graphical construction, however carefully reinvestigated, did not 
permit of our considering the curve to be anything but a straight 
line. . . Even if it is not absolutely true, it exemplifies the extraordinary 
power of such integrals of J products [that is, (5)] to give extremely 
close approximations to such simple forms as horizontal lines." 

This conjecture was investigated in detail in [Fet63] wherein the nonlinearity was 
ffist rigorously established. This work and various more recent papers highlight the 
difficulty of computing the underlying Bessel integrals. 

Remark 8. Recall from Example 4 that the asymptotic behaviour of Pn at zero is 
determined by the poles of the moments Wn{s). To obtain information about the 
behaviour of p„(x) as x — )■ n^, we consider the "reversed" densities p„(x) = pn{n — x) 
and their moments Wn{s). For non-negative integers k, 

Wn{k) = ^ X'pnix) dx = J^{n- x)Vn(x) = 5Z Q 1)'^"'' (j ) 
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On the other hand, we can find a recurrence satisfied by the Wn{s) as follows: a 
differential equation for the densities Pn{x) is obtained from Theorem 2 by a change 
of variables. The Mellin transform method as described in Example 1 then provides a 
recurrence for the moments Wn{s). We next apply the same reasoning as in [BSWll] 
to obtain information about the pole structure of Wn{s). It should be emphasized 
that this involves knowledge about initial conditions in term of explicit values of 
initial moments VF„(2/c). 

For instance, in the case n = 4, we find that the moments H/4(s) have simple poles 
at — |,— |,— which predicts an expansion of p^i^x) as given in Theorem 5. 

For n = 5, we learn that W^ls) has simple poles at s = —2, —3, —4, .... It then 
follows, as for (52), that P5{x) = J^kLo^^,!: ~ 5)'^^^ for x ^ 5 and close to 5. The 
fs^fc are the residues of ^^5(5) at s = —k — 2. 



6. Derivative evaluations of Wn 

As illustrated by Theorem 6, the residues of Wn{s) are very important for studying 
the densities p„ as they directly translate into behaviour of p„ at 0. The residues 
may be obtained as a linear combination of the values of Wn{s) and VF^(s). 

Example 5 (Residues of Wn)- Using the functional equation for W^ls) and L'Hopital's 
rule we find that the residue at s = — 2 can be expressed as 

(54) Re=.,(,y3) - + 

y 

This is a general principle and we likewise obtain for instance: 

16 + 1140W^^(0) - 8041V^(2) + 64iy^(4) 



(55) Res_2(iy5 

(56) Res_4(Vr5 



225 

26Res_2(Pr5) - 16 - 20W^(0) + 4iy^(2) 



225 

In the presence of double poles, as for W4, 



(57) lim {s + 2yW4{s) 
and for the residue: 

(58) Res_,(»y = 8+18^n(0)-3W--(2)+4tn'(0)-H?(2)^ 
Equations (57, 58) are used in Example 4 and each unknown is evaluated below. 



We are therefore interested in evaluations of the derivatives of Wn for even argu- 
ments. 
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Example 6 (Derivatives of W3 and W4). Differentiating the double integral for 
1^3(5) (7) under the integral sign, we have 

H^^(O) = ^ /" f log(4 sm{7iy) cos(27rx) + 3-2 cos(27r?/)) dx dy. 
2 Jo Jo 

Then, using 

J log(a + 6cos(27ra;)) dx = log (^a + Va^ — for a > 6 > 0, 



we deduce 



(59) W'M = / log(2sin(7r|/)) dy = - CI 

Jl/6 V 

where CI denotes the Clausen function. Knowing as we do that the residue at s = —2 
is 2/(v^vr), we can thus also obtain from (54) that 

'-3P)^2.ia(^)-^. 

In like fashion, 

^4(0) — — 9 / / log (3 + 2 cos X + 2 cos y + 2 cos{x — y)) dx dy 
8^ Jo Jo 

The final equality will be shown in Example 8. Note that we may also write 

^3(0) = — 9 / / log(3 + 2 cosx + 2 cosy + 2 cos(x — y)) dx dy. 
Stt^ Jo Jo 

The similarity between ^^3(0) and W^4(0) is not coincidental, but comes from applying 

log ((a + cos27rx) + (6 + sm27rx) ) dx = | ^ ^ ^^^^^^.^^ 

to the triple integral of 14^4(0). As this reduction breaks the symmetry, we cannot 
apply it to 1^5(0) to get a similar integral. 

In general, differentiating the Bessel integral expression 
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obtained by David Broadhurst [Bro09] and discussed in [BSWll], under the integral 
sign gives 

2 » » -ra-l, 



W'S^) = n log - - 7 J^-\x)Ux) dx 



X . 

oo 



/>oo 

(62) =log(2)-7-n/ \og{x)J^-^{x)Ji{x)dx, 

Jo 

where 7 is the Euler- Mascheroni constant, and 

C(0) =nj^^ (^log 0^ - 7) Jr\^)M^) dx. 

POO 

= (log(2) - l)W4-l) - / log(x) Jo"(x) dx, 

Jo 

= / - Jr'(x) Ji(x) (1 - 7 - log(2a;)) dx. 
Jo ^ 

Remark 9. We may therefore obtain many identities by comparing the above equa- 
tions to known values. For instance, 



Likewise 



and 



3 / log{x)J^{x)Ji{x) dx = log(2) - 7 - - CI ( 

Jo TT V 



37- 



Example 7 (Derivatives of W5). In the case n = 5, 

ly^(O) = 5 ^ (^og 0^ - 7^ Jo (t) Ji(t) dt ^ 0.54441256 



with similar but more elaborate formulae for Wl^{2) and ^^5(4). Observe that in 
general we also have 

dx 

X 

which is useful numerically. 



(63) 



iy;(0) = log(2) - 7 - - 1) - - r J^ix) 

Jo ^ Ji 



In fact, the hypergeometric representation of and PV4 first obtained in [Cra09] 
and recalled below also makes derivation of the derivatives of W3 and possible. 

Corollary 2 (Hypergeometric forms). For s not an odd integer, we have 
(64) 



W^(s) = —r— tan 



1 /TTSX / S^' ^1^1 



2' 2' 2 



s s 



2' 2' 2 
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and, if also Res > —2, we have 
(65) 



WJs] 



1 /ns 
—— tan — 
V 2 



2 



4-^3 



1 1 1 £ + 1 

2 ' 2 ' 2 ' 2 ^ 
s+3 s+3 s+3 

2 ' 2 ' 2 



4-^3 



2' 2' 2' 



1,1, 



s-1 
2 



Example 8 (Evaluation of ^^^^(O) and Wi{0)). If we write (64) or (65) as W^is) = 
/i(s)Fi(s) + /2(s)F2(s), where Fi, F2 are the corresponding hypergeometric functions, 
then it can be readily verified that /i(0) = j^i^) = -^2(0) — 0- Thus, differentiating 
(64) by appealing to the product rule we get: 



TT 



"3-^2 



111 
2' 2' 2 
3 3 
2' 2 



71 \ 3 



The last equality follows from setting = 71/6 in the identity 



2sin(^)3F2 



111 

2' 2' 2 
3 3 
2' 2 



sin^ 9) = CI (2 0) + 2 log (2 sin ( 



(66) 



Likewise, differentiating (65) gives 

4 „ 



71^ 



74-^3 



1 1 1 - 

2 ' 2 ' 2 ' " 
3 3 3 
2' 2' 2 



703) 
27r2 ^ 



thus verifying (60). In this case the hypergeometric evaluation 

1 1 1 



4-^3 



2' 2' 2' ■ 
3 3 3 
2' 2' 2 



n=0 



7 



(2n + l)3 8 



C(3), 



is elementary. 

Differentiating (64) at s = 2 leads to the evaluation 



3-^2 



/I 1 1 

1 2' 2' 2 





27 / 


1 5 5 

^ 2' 2 







CI 



71 



while from (65) at s = 2 we obtain 
(67) wm 



14C(3) - 12 



71^ 







Thus we have enough information to evaluate (57) (with the answer 3/(271^)). 

Note that with two such starting values, all derivatives of ^^3(5) or W^i^s) at even 
s may be computed recursively. 
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We also note here that the same technique yields 

12 TT ^ 42" (2n + 1)2 

n=0 ^ ^ 



(68) ^3(o)-^;-!f:S) 



^ ' 12^ TT \3) vr ^ 42" (2n + 1)2 

n=0 ^ ^ 



and, quite remarkably, 

(^0) M/aO)-^2+ ^2 +vr2 2- (2n + l)3 

n=0 ^ ^ 

_ 24Li4 (I) - 18C(4) + 21C(3) log(2) - 6C(2) log^(2) + log^(2) 

7r2 

where the very final evaluation is obtained from results in [BZB08, §5]. Here Li4(l/2) 
is the polylogarithm of order 4, while Hn := 7 + \E'(n + 1) denotes the nth harmonic 
number, where \E' is the digamma function. So for non-negative integers n, we have 
explicitly Hn = '^k=i ^/^^ before, and 



^n+i/2 = 2^^^-21og(2). 

k=l 



An evaluation of W^'{0) in terms of polylogarithmic constants is given in [BSll]. 
In particular, this gives an evaluation of the sum on the right-hand side of (68). 

Finally, the corresponding sum for W^{2) may be split into a telescoping part and 
a part involving W^{0). Therefore, it can also be written as a linear combination 
of the constants used in (70). In summary, we have all the pieces to evaluate (58), 
obtaining the answer 9 log(2)/(27r2). 

6.1. Relations with Mahler measure. For a (Laurent) polynomial f{xi, . . . ,x„), 
its logarithmic Mahler measure, see for instance [RVTV04], is defined as 









• / log 


/o 


Jo 



m{f)= ... log f(e'^''\...,e'^''-) dh---dtn. 



Recall that the sth moments of an n-step random walk are given by 



WJs) 



n 



k=l 



dti - ■ ■ dtn = \\Xi + . . . + Xn\ 



'0 JO 

where || ■ \\p denotes the p-norm over the unit n-torus, and hence 

W!^{0) = m{xi + . . . + Xn) = m{l + xi + . . . + x„_i). 
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Thus the derivative evaluations in the previous sections are also Mahler measure 
evaluations. In particular, we rediscovered 

ly^(O) = ^ Cl (^) = L'(x-3, -1) = m(l + xi + X2), 

along with 



w'M = 4V = + X1 + X2 + X3) 



703) 

which are both due to C. Smyth [RVTV04, (1.1) and (1.2)] with proofs first published 
in [BoySl, Appendix 1]. 

With this connection reahzed, we find the following conjectural expressions put 
forth by Rodriguez- Villegas, mentioned in different form in [Fin05] , 



V / 15 ^ 5/2 ,00 



(71) W'M=[-^,] I {r/^(e->^(e-^*)+r^^(e-V(e-^'*)}t'dt 



and 



7 



3 



00 

2/ -t\^2t -2t\^2t -U\^2l -&t\ J.4 



(72) W',{Q) = ( ^ ) / r^\e-'W{e-'''W{e--''W{e-'') dt, 







where 77 was defined in (42). We have confirmed numerically that the evaluation of 
14^5(0) in (71) holds to 600 places. Likewise, we have confirmed that (72) holds to 80 
places. Details of these somewhat arduous confirmations are given in [BBIO]. 

Differentiating the series expansion for Wn{s) obtained in [BNSW09] term by term, 
we obtain 

m=l k=0 ^ ' 

On the other hand, from [RVTV04] we find the strikingly similar 

ra=2 k=0 ^ ^ 

Finally, we note that Wni^s) itself is a special case of zeta Mahler measure as 
introduced recently in [Aka09]. 

7. New results on the moments Wn 
From [BBBG08] equation (23), we have for A; > even, 

oA:+3/2 fco 
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where Io(t),Ko(t) denote the modified Bessel functions of the first and second kind, 
respectively. 

Similarly, [BBBG08] equation (55) states that for A; > even, 

Ak+2 poo 

Equation (75) can be formally reduced to a closed form as a 3F2 (for instance by 
Mathematica) . At k = ±1, the closed form agrees with VF3(±1). As both sides of 
(75) satisfy the same recursion ([BBBG08] equation (8)), we see that it in fact holds 
for all integers k > —2. 

In the following we shall use Carlson's theorem ([Tit39]) which states: 
Let f be analytic in the right half-plane Kez ^ and of exponential type with the 
additional requirement that 

\f{z)\ ^ Me'^l^l 

for some d < tt on the imaginary axis Re 2; = 0. // f{k) = for k = 0,1,2, .. . then 
f{z) = identically. We then have the following: 

Lemma 1. Equation (75) holds for all k with Kek > —2. 

Proof. Both sides of (75) are of exponential type and agree when = 0, 1, 2, . . .. The 
standard estimate shows that the right-hand side grows like el?^!''/^ on the imaginary 
axis. Therefore the conditions of Carlson's theorem are satisfied and the identity 
holds whenever the right-hand side converges. □ 

Using the closed form given by the computer algebra system, we thus have: 

Theorem 11 (Single hypergeometric for W^ls)). For s not a negative integer < — 1, 

4 



os+3/2 Yd I ^/o] 



2 / s+2 3+2 s+2 

2 ' 2 ' 2 

1 £+3 
2 



Turning our attention to negative integers, we have for ^ an integer: 
(78) w,(-2k-l)=-,(^^) I t^^-Mtfit, 

because the two sides satisfy the same recursion ([BBBG08, (8)]), and agree when 
k = 0,l ([BBBG08, (47) and (48)]). 

Remark 10. Equation (78) however does not hold when k is not an integer. Also, 
combining (78) and (75) for Wsi^—l), we deduce 

/oo o roo 2 poo 

Ko{tfmdt = -j^j^ Ko{tfdt = ^j^ utfdt. 
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From (78), we experimentally determined a single hypergeometric for ^^3(5) at 
negative odd integers: 



Lemma 2. For k ^ an integer, 



W.(-2k - 11 



3-^2 



111 
2' 2' 2 

k + l,k + l 



Proof. It is easy to check that both sides agree at k = 0, 1. Therefore we need 
only to show that they satisfy the same recursion. The recursion for the left-hand 
side implies a contiguous relation for the right-hand side, which can be verified by 
extracting the summand and applying Gosper's algorithm ([PWZ06]). □ 

The integral in (78) shows that Ws{—2k — 1) decays to rapidly - very roughly 
like as — )■ CX3 - and so is never when k is an integer. 

To show that (76) holds for more general k required more work. Using Nicholson's 
integral representation in [Wat41], 

»7r/2 



the integral in (76) simplifies to 



Ko{2ts'm a) da. 



(79) 



2 



/'7r/2 poo 

/ / t''-^^Ko{tyKo{2tsma)dtda. 
Jo Jo 



The inner integral in (79) simplifies in terms of a Meijer G-function; Mathematica 
is able to produce 



-y3,2 



I sm 



k+2 ^ ^3,3 



1 _1 1 
'2' 2' 2 

0,0,0 



sin^ a 



which transforms to 



1 sm 



7^ ^2,3/1,1,1 

k+2 „ *-^3,3 (331 



a 



2' 2' 2 



sin^ a 



Let t = sin^ a in the above, so the outer integral in (79) transforms to 



(80) 



'IT 

16 



r^(i-t)-G^ 





1,1 




(3' 


3 1 




V 2' 


2' 2 





We can resolve this integral by applying the Euler-type integral 



'1) 



r"(i -t) 



a— 6—1 rimm 



G;:: l^t) dt = r(a - h) G^X,\i 



a, c 

h 



Indeed, when k = —1, the application of (81) recovers the Meijer G representation 
of Wii-l) ([BSWll]); that is, (76) holds for k = -1. 
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When k = 1, the resulting Meijer G-function is 

x-y2,4 [ 2, 1, 1, 1 

4,4 1 3 3 1 3 
V 2' 2' 2' 2 

to which we apply Nesterenko's theorem ([Nes03]), deducing a triple integral (up to 
a constant factor) for it: 



1 rl 



^ ^ dxdydz. 



lo Jo Jo ^ y{l-y){l-z){l-x{l-yz)y^ 
We can reduce the triple integral to a single integral, 

"1 8E'{t){{l + e)K'{t)-2E'{t)) 



dt. 



Now applying the change of variable t h-?- (1 — t)/(l + t), followed by quadratic 
transformations for K and E, we finally get 



which is, indeed, (a correct constant multiple times) the expression for ^4(1) which 
follows from Section 3.1 in [BSWll]. 

We finally observe that both sides of (76) satisfy the same recursion ([BBBG08] 
equation (9)), hence they agree for k = 0, 1,2, . . .. Carlson's theorem applies with 
the same growth on the imaginary axis as in (75) and we have proven the following: 

Lemma 3. Equation (76) holds for all k with Kek > —2. 

Theorem 12 (Alternative Meijer G representation for W^ls)). For all s, 



n2s+l / 111 s+3 

2 Ff'l /'o _L 0^^2 4,4 I s+2 s+2 s+2 1 



^2) W^4(s) = ^^-^7 ^G^'^' ' ' ' 2 

^ ^ ^ vr2r(l(. + 2))2 • V 2 ' 2 ' 2 '2 



Proof. Apply (81) to (80) for general k, and equality holds by Lemma 3. □ 

Note that Lemma 3 combined with the known formula for W^l—l) in [BSWll] 
gives 

4 



y]-3 



Koityioit)dt= / Jo(t)Mt 



Armed with the knowledge of Lemma 3, we may now resolve a very special but 
central case (corresponding to n = 2) of Conjecture 1 in [BSWll]. 
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Theorem 13. For integer s, 

(83) W,{s) = Y,r- W,{s-2j). 

j=o V J / 

Proof. In [BNSW09] it is shown that both sides satisfy the same three term recur- 
rence, and agree when s is even. Therefore, we only need to show that the identity 
holds for two consecutive odd values of s. 
For s = — 1, the right-hand side of (83) is 

°° /_i /oX 2 oo r)2~2j poo 

j=o \ J / j=o J- JO 

upon using (78), and after interchanging summation and integration (which is justi- 
fied as all terms are positive), this reduces to 



TT' 



4 rco 
~3 



Ko{tflo{t) dt, 



which is the value for W^l—l) by Lemma 3. 

We note that the recursion for ^4(5) gives the pleasing reflection property 



W4i-2k - 1) 2^'= = 1^4(2^ - 1). 

In particular, W4(— 3) = ^^^4(1). Now computing the right-hand side of (83) at 
s = —3, and interchanging summation and integration as before, we obtain 

(^~jyW,{-3-2j) = ± t'Ko{tflo{t)dt = ^W,{1) = W,{-3). 
Therefore (83) holds when s = — 1, —3, and thus holds for all integer s. □ 
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Appendix. A family of combinatorial identities 



DON ZAGIERi 

The "collateral result" of Djakov and Mityagin, [DM04, DM07], is the pair of 
identities 

E ^I(m^-^D=^.(1^3^...,(2m-ln, 

— m<ii<---<ife<m s=l 
*2-«iv,«fc-*fe-i>2 

k 

J2 Hi"^ - ^«)(^ + is-l)= (Tk{2^, 4', . . . , (2m - 2f) , 

l—m<ii<---<if;<m s=l 
12-41, .■■,«fc-*fe-i>2 

where m and k are integers with m ^ k ^ and ak denotes the fcth elementary 
sjTiimetric function. By setting js = ig + m in the first sum and = + m — 1 in 
the second, we can rewrite these formulas more uniformly as^ 



ak{l\ 32, . . . , (M - 1)2) if M is even, 
crfc(22,42,...,(M- 1)2) if Mis odd. 



(84) FmAM) ^ 

where FM,k{X) is the polynomial in X (non-zero only if M ^ 2A; ^ 0) defined by 



(85) FM,k{X)= E n^'^(^-^'^)- 

0<ii<-<jfe<Af s=l 

i2-ii,...,jfc-jfc-i>2 

The advantage of introducing the free variable X in (85) is that the functions FM.k{X) 
satisfy the recursion 

(86) FM+i,k+i{X) - FM,k+i{X) = M{X-M) Fm-iA^), 

because the only paths that are counted on the left are those with < ji < ■ ■ ■ < 

jk < jk+i = M. 

It is also advantageous to introduce the polynomial generating function 



$M = $m(X,m)= E (-l)'^M,fc(X) 

0s;fcs£M/2 



^The original note is unchanged. 
^Note that (84) is precisely Theorem 4. 



DENSITIES OF SHORT UNIFORM RANDOM WALKS 29 

the first examples being 

$0 = 1, $i = M, $2 = ?i^-(X-l), $3 = ^3- (3X-5)m, 
$4 = - (6X - 14)^2 + (3X2 _ + 9) , 
$5 = - (lOX - 30)m^ + (ISX^ - SOX + 89) , 

$g = _ (15X - 55)m^ + (45X2 - 300X + 439)^2 _ (isx^ - 135X2 + 345X - 225) . 
In terms of tliis generating function, the recursion (86) becomes 
(87) <l>M+i = u - M(X - M) ^M^i 

and the identity (84) to be proved can be written succinctly as 



^m{M,u)= W {u-X). 



\\\<M 
X^M (mod 2) 

Denote by Pm{u) the polynomial on the right-hand side of (88). Looking for other 
pairs (M, X) where $m(-^;W) has many integer roots, we find experimentally that 
this happens whenever M — X is a non-negative integer, and studying the data more 
closely we are to conjecture the two formulas 

(89) $m(M -n,u) = -Y,[-)PM{u-n + 2j) (M, n ^ 0) 

j=0 

(a generalization of (88)) and 

(90) $iv/+n(M, u) = $m(M, u) $„(-M, u) (M, n^O). 

Formula (90) is easy to prove, since it holds for n = trivially and for n = 1 by (87) 
and since both sides satisfy the recursion = uyn+n{M +n) Vn-i for n = 1, 2, . . . 
by (87). On the other hand, combining (88), (89) and (90) leads to the conjectural 
formula 



or, renaming the variables, 

(91) M »m(- + y + 1, y - = (;) (,/_ ^ ■ 

To prove this, we see by (87) that, denoting by Gm = GM{x,y) the expression on 
the right, it suffices to prove the recursion (M -|- 1)Gm+i = {y — x)Gm + {M — x — 
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y — 1)Gm-i- This is an easy binomial coefficient identity, but once again it is easier 
to work with generating functions: the sum 

oo 

(92) g{x, y- T) := Gm(x, y) = (1 - Ty (1 + T)^ 

A/=0 

satisfies the differential equation 

1 dQ y X 
gdf ~ 1 + T ~ 1 - T 

or 

§ = (y-x)^+(rA_^_^^ ^' 

and this is equivalent to the desired recursion. 

We can now complete the proof of (84). Rewriting (92) in the form 

we find that, for 1 ^ j ^ M, 

^ <I>m(M, M + 1 - 2j) = coeff^M ((1 - Ty-' (1 + T)^^--') = 

and hence that the polynomial on the left-hand side of (88) is divisible by the poly- 
nomial on the right, which completes the proof since both are monic of degree M 
in u. 
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